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The Commissioner of Patents and Trademarks 
Washington, D.C. 20231 



PRELIMINARY AMENDMENT 



Dear Sir: 



By way of preliminary amendment, please amend the above identified application as set 
forth below. 

In the Specification; 

Please cancel pages 1-3 of the specification, in their entirety, in favor of a clean form 
of pages 1-3 of the specification, without any markings thereon, accompanying this response. 
Also accompanying this response is a copy of the original pages of the specification which show 
the addition(s) (in red ink) and the deletion(s) (by strikeout) to the canceled specification pages. 
The amended specification pages do not include any new subject matter and only include the 
changes which are indicated on the marked-up copy of the specification pages. Please enter 
the replacement specification pages into the record of this case. 

In the Claims : 

Please cancel original claims 1 - 16, as well as any Chapter II amended 
claims, in favor of new claims 17 -33. 



17. (NEW) A method of imaging in which a de-convolution process is applied to 
the image-domain results of the object-scan to derive therefrom the respective point- or line- 
spread function of at least one object-discontinuity, and to derive from said function the 
location in the image domain of the respective discontinuity, wherein the de-convolution 
process is carried out using sub-pixel sampling. 

18. (NEW) The method according to claim 17, wherein the location of the 
respective discontinuity is derived from the mid-point of the full-width half-maximum of the 
function. 

19. (NEW) The method according to claim 17, wherein the function is correlated 
with the image-domain results of the transfer for enhancement of spatial resolution of the 
imaging of the respective discontinuity. 

20. (NEW) The method according to claim 19, wherein the enhancement of 
spatial resolution of the imaging of the respective discontinuity involves transfer of sub-pixels 
within the image-domain results of the respective discontinuity, the sub-pixels being 
transferred within their respective image-domain results from one side to the other of the 
location for edge-image definition. 

21 . (NEW) The method according to claim 17, wherein the de-convolution 
process is carried out using least-squares running filtering. 

22. (NEW) The method according to claim 17, wherein an edge-contour of the 
object is defined in the image domain using de-convolution processing as aforesaid. 

23. (NEW) The method according to claim 22, wherein the area of the object- 
image within the edge-contour is determined. 

24. (NEW) The method according to claim 22, wherein the volume of the object- 
image within the edge-contour is determined. 
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25. (NEW) The method according to claim 22, wherein the object-scan is a 
magnetic resonance (MR) scan, values of relaxation times T,and T 2 are derived for the 
object-image within said contour, and these values are used to identify from stored data, 
types of material involved in the scanned object. 

26. (NEW) The method according to claim 25, wherein density values for the 
identified material types are derived from further stored data. 

27. (NEW) The method of imaging according to claim 17, wherein corresponding 
computed tomography (CT) and magnetic resonance (MR) scans of the same part of an 
object are derived, the scans are related to one another for correlation of one to the other 
positionally with respect to said pari using the de-convolution process, and imaging of said 
part of the object is provided in accordance with the MR scan as modified spatially in 
dependence upon the CT contrast numbers applicable to the corresponding, correlated 
positions of the CT scan. 

28. (NEW) The method according to claim 27, wherein geometric correction is 
applied to the imaging derived from the MR scan, in accordance with stored data. 

29. (NEW) An imaging system comprising mechanism for performing a de- 
convolution process on the image-domain results of an object-scan to derive therefrom the 
respective point- or line-spread function of at least one object-discontinuity, and mechanism 
to derive from the function the location in the image domain of the respective discontinuity, 
wherein the de-convolution process is carried out using sub-pixel sampling. 

30. (NEW) The system according to claim 29, wherein the location of the 
respective discontinuity is derived from the mid-point of the full-width half-maximum of the 
function. 

31 . (NEW) The system according to claim 29, wherein the function is correlated 
with the image-domain results of the transfer for enhancement of spatial resolution of the 
imaging of the respective discontinuity. 
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32. (NEW) The system according to claim 31 , wherein the enhancement of spatial 
resolution of the imaging of the respective discontinuity involves transfer of sub-pixels within 
the image-domain results of the respective discontinuity, the sub-pixels being transferred 
within their respective image-domain results from one side to the other of the location for 
edge-image definition. 

33. (NEW) The system according to claim 29, wherein the de-convolution process 
is carried out using least-squares running filtering. 
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REMARKS 

Accompanying this response, please find replacement pages and marked-up pages 
of the specification which overcome some informalities noted in the specification on file. The 
undersigned avers that the enclosed replacement page(s) of the specification do not contain 
any new matter. 

In the event that there are any fee deficiencies or additional fees are payable, please 
charge the same or credit any overpayment to our Deposit Account (Account No. 04-0213). 

Respectfully submitted, 




Michael J. Bujold^ 
Customer No. 
Davis & Bujold, P.L.L.C^ 
Fourth Floor 

500 North Commercial Street 
Manchester NH 03101-1151 
Telephone 603-624-9220 
Facsimile 603-624-9229 
E-mail: patent@davisandbujold.com 
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Imaging 

This invention relates to imaging and in particular to 
methods and systems for image enhancement. 

Imaging involves transfer from the object domain into the 
image domain, but owing to limiting factors such as the 
finite size of energy source, detector size, sampling 
frequency, display density, software filter function, and 
possibly partial-volume effects experienced with some 
imagers, an infinitely fine delta function in the object 
domain cannot be faithfully reproduced in the image 
domain ♦ Instead, a smeared-out image, or point-spread 
function (PSF) , is observed. Similarly, an infinitely 
sharp edge-response function (ERF) in the object domain 
becomes a smeared-out ERF in the image domain. The 
smearing effect becomes more intense as the adjacent ERFs 
of discontinuities or contrast profiles get closer to 
each other. 

It is an object of the present invention to provide a 
method and system by which the above problem can be at 
least .partly overcome. 

SOMm^ of nth /WeMtW 

According to one aspect of the present invention there is 
provided a method wherein a de-convolution process is 
applied to the image-domain results of an object-scan to 
derive therefrom the respective point- or line-spread 
function effective in the object- to image-domain 
transfer of one or more object-discontinuities, and to 
derive from said function the location in the image 
domain of the respective discontinuity. . 

According to another aspect of the invention there is 
provided an imaging system comprising means for 
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performing a de-convolution process on the image-domain 
results of an object-scan to derive therefrom the 
respective point- or line-spread function effective in 
the object- to image-domain transfer of one or more 
5 object-discontinuities, and means to derive from said 
function the location in the image domain of the 
respective discontinuity. 

The method and system of the invention enable the 
10 location of the respective discontinuity in the image 

domain, to be established with a high degree of accuracy. 
This is critical to image definition free of any 
substantial smearing, and to this end the location of the 
respective discontinuity may be derived to sub-pixel 
15 accuracy simply from the mid-point of the full-width 
half -maximum of said function. 

The said function may be correlated with the image-domain 
results of said transfer for enhancement of spatial 

20 resolution of the imaging of the one or more 

discontinuities. This enhancement may involve transfer 
of sub-pixels within the image-domain results of the 
respective one or more discontinuities, the sub-pixels 
being transferred within their respective image-domain 

25 results from one side to the other of said location. 

The de-convolution process may be carried out using 
least-squares running filtering. 

filler b£ScgipT7ou of the h&m^Ss 

30 An imaging method and system according to the present 

invention will now be described, by way of example, with 
reference to the accompanying drawings, in which: 

Figure 1 illustrates schematically the method and system 
35 of the invention; 
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Figure 2 illustrates features of processing performed in 
the method and system of Figure 1; 

Figure 3 shows results achieved from use of the method 
5 and system of Figure 1; 

Figure 4 shows to an enlarged scale a section of the 
contour of an image profile depicted in Figure 2; 

10 Figures 5 and 6 are a plan view and sectional end- 
elevation of a couch-top used in the method and system of 
Figure 1; and 

Figure 7 provides illustrates of a convolution operation, 
15 as a basis for a mathematical model of de-convolution 
processing in accordance with the method and system of 
Figure 1. 

The method and system to be described with reference to 
20 Figure 1 utilise MR scanning for medical diagnostic and 
treatment-planning purposes. In principle and in the 
general techniques described, the method and system of 
the invention can be used in other applications of MR 
scanning and also in circumstances where other scanning 
25 techniques are utilised* Furthermore, although both 
structure and function are represented by discrete 
x boxes' 1 to 19 in Figure 1, the method and system are to 
a substantial extent manifest in programmed digital data- 
processing operations. 

30 

Referring to Figure 1, data derived in accordance with 
conventional operation of an MR scanner 1 is processed 
for imaging purposes within a processor 2 . The output of 
the processor 2 is used to provide a display 3 , and from 
35 this is subject to post-imaging processing 4. The post- 
imaging processing 4 includes the facility for selecting 
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Imaging 

FIELD OF THE INVENTION 

5 This invention relates to imaging and in particular to 
methods and systems for image enhancement. 

BACKGROUND OF THE INVENTION 

Imaging involves transfer from the object domain into the 
image domain, but owing to limiting factors such as the 

10 finite size of energy source, detector size, sampling 

frequency, display density, software filter function, and 
possibly partial -volume effects experienced with some 
imagers, an infinitely fine delta function in the object 
domain cannot be faithfully reproduced in the image 

15 domain. Instead, a smeared-out image, or point-spread 
function (PSF) , is observed. Similarly, an infinitely 
sharp edge-response function (ERF) in the object domain 
becomes a smeared-out ERF in the image domain. The 
smearing effect becomes more intense as the adjacent ERFs 

20 of discontinuities or contrast profiles get closer to 
each other. 

It is an object of the present invention to provide a 
method and system by which the above problem can be at 
25 least. .partly overcome, 

SUMMARY OF THE INVENTION 

According to one aspect of the present invention there is 
provided a method wherein a de-convolution process is 
applied to the image-domain results of an object-scan to 
30 derive therefrom the respective point- or line-spread 
function effective in the object- to image-domain 
transfer of one or more object-discontinuities, and to 
derive from said function the location in the image 
domain of the respective discontinuity. _ 

35 

According to another aspect of the invention there is 
provided an imaging system comprising means for 
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performing a de-convolution process on the image-domain 
results of an object-scan to derive therefrom the 
respective point- or line-spread function effective in 
the object*- to image-domain transfer of one or more 
5 object-discontinuities, and means to derive from said 
function the location in the image domain of the 
respective discontinuity. 

The method and system of the invention enable the 
10 location of the respective discontinuity in the image 

domain, to be established with a high degree of accuracy. 
This is critical to image definition free of any 
substantial smearing r and to this end the location of the 
respective discontinuity may be derived to sub-pixel 
15 accuracy simply from the mid-point of the full-width 
half -maximum of said function. 

The said function may be correlated with the image-domain 
results of said transfer for enhancement of spatial 

20 resolution of the imaging of the one or more 

discontinuities. This enhancement may involve transfer 
of sub-pixels within the image-domain results of the 
respective one or more discontinuities, the sub-pixels 
being transferred within their respective image-domain 

25 results from one side to the other of said location. 

The de-convolution process may be carried out using 
least-squares running filtering. 
BRIEF DESCRIPTION OF THE DRAWINGS 
30 An imaging method and system according to the present 

invention will now be described, by way of example, with 
reference to the accompanying drawings, in which: 

Figure l illustrates schematically the method and system 
35 of the invention; 



WO 0022573 



PCT/GB99/03417 



3 . 

Figure 2 illustrates features of processing performed in 
the method and system of Figure 1; 

Figure 3 shows results achieved from use of the method 
5 and system of Figure 1; 

Figure 4 shows to an enlarged scale a section of the 
contour of an image profile depicted in Figure 2; 

10 Figures 5 and 6 are a plan view and sectional end- 
elevation of a couch-top used in the method and system of 
Figure 1; and 

Figure 7 provides illustrates of a convolution operation, 
15 as a basis for a mathematical model of de-convolution 
processing in accordance with the method and system of 
Figure 1. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The method and system to be described with reference to 

20 Figure 1 utilise MR scanning for medical diagnostic and 
treatment-planning purposes. In principle and in the 
general techniques described, the method and system of 
the invention can be used in other applications of MR 
scanning and also in circumstances where other scanning 

25 techniques are utilised. Furthermore, although both 
structure and function are represented by discrete 
1 boxes' 1 to 19 in Figure 1, the method and system are to 
a substantial extent manifest in programmed digital data- 
processing operations. 

30 

Referring to Figure 1, data derived in accordance with 
conventional operation of an MR scanner 1 is processed 
for imaging purposes within a processor 2 . The output of 
the processor 2 is used to provide a display 3, and from 
35 this is subject to post-imaging processing 4. The post- 
imaging processing 4 includes the facility for selecting 
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This invention relates to imaging and in particular to 
5 methods and systems for image enhancement. 

Imaging involves transfer from the object domain into the 
image domain , but owing to limiting factors such as the 
finite size of energy source, detector size, sampling 

10 frequency, display density, software filter function, and 
possibly partial-volume effects experienced with some 
imagers, an infinitely fine delta function in the object 
domain cannot be faithfully reproduced in the image 
domain. Instead, a smear ed-out image, or point-spread 

15 function (PSF) , is observed. Similarly, an infinitely 
sharp edge-response function (ERF) in the object domain 
becomes a smear ed-out ERF in the image domain. The 
smearing effect becomes more intense as the adjacent ERFs 
of discontinuities or contrast profiles get closer to 

2 0 each other. 

It is an object of the present invention to provide a 
method and system by which the above problem can be at 
least partly overcome. 

25 

According to one aspect of the present invention "there is 
provided a method in which a de-convolution process is 
applied to the image-domain results of an object-scan to 
derive therefrom the respective point- or line-spread 

3 0 function of one or more object-discontinuities, and to 

derive from said . function the location in the image 
domain of the respective discontinuity, wherein the de- 
convolution process is carried out using sub-pixel 
sampling. 



According to another aspect of the invention there is 
provided an imaging system comprising means for 
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performing a de-convolution process on the image-domain 
results of an object-scan to derive therefrom the 
respective point- or line-spread function of one or more 
object-discontinuities, and means to derive from said 
function the location in the image domain of the 
respective discontinuity , wherein the de-convolution 
process is carried out using sub-pixel sampling. 

The method and system of the invention enable the 
location of the respective discontinuity in the image 
domain, to be established with a high degree of accuracy. 
This is critical to image definition free of any 
substantial smearing, and to this end the location of the 
respective discontinuity may be derived to sub-pixel 
accuracy simply from the mid-point of the full-width 
half -maximum of said function. 

The said function may be correlated with the image-domain 
results of said transfer for enhancement of spatial 
resolution of the imaging of the one or more 
discontinuities. This enhancement may involve transfer 
of sub-pixels within the image-domain results of the 
respective one or more discontinuities, the sub-pixels 
being transferred within their respective image-domain 
results from one side to the other of said location for 
edge-image definition. 

The de- convolution process may be carried out using 
least-squares running filtering. 

An imaging method and system according to the present 
invention will now be described, by way of example, with 
reference to the accompanying drawings, in which: 

Figure 1 illustrates schematically the method and system 
of the invention; 
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Figure 2 illustrates features of processing performed in 
the method and system of Figure 1; 

Figure 3 shows results achieved from use of the method 
5 and system of Figure 1; 

Figure 4 shows to an enlarged scale a section of the 
contour of an image profile depicted in Figure 2; 

10 Figures 5 and 6 are a plan view and sectional end- 
elevation of a couch-top used in the method and system of 
Figure 1; and 

Figure 7 provides illustrates of a convolution operation, 
15 as a basis for a mathematical model of de-convolution 
processing in accordance with the method and system of 
Figure l. 

The method and system to be described with reference to 
20 Figure 1 utilise MR scanning for medical diagnostic and 
treatment-planning purposes. In principle and in the 
general techniques described, the method and system of 
the invention can be used in other applications of MR 
scanning and also in circumstances where other scanning 
25 techniques are utilised ♦ Furthermore, although both 
structure and function are represented by discrete 
11 boxes' 1 to 19 in Figure 1, the method and system are to 
a substantial extent manifest in programmed digital data- 
processing operations. 

30 

Referring to Figure 1, data derived in accordance with 
conventional operation of an MR scanner 1 is processed 
for imaging purposes within a processor 2 . The output of 
the processor 2 is used to provide a display 3 , and from 
35 this is subject to post-imaging processing 4. The post- 
imaging processing 4 includes the facility for selecting 
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a region of the display 3 for more-detailed and closer 
inspection. 

To the extent the imaging method and system of Figure 1 
have so far been described, they are conventional, and it 
is in further processing 5 of the image data of the 
selected region of interest obtained by the post-imaging 
processing 4 that a step forward from what is already 
known is achieved* More particularly, the further 
processing 5 is operative to define more clearly the true 
edges or boundaries of image contour (s) in the selected 
region of interest, and to enhance the accuracy of the 
imaging of those contours- 

The definition and accuracy of transfer of features from 
the object domain scanned by the scanner 1, to the image 
domain manifest in the post-imaging processing 4, is 
limited by many factors. The limitations arise from 
within the scanner 1 itself (in particular the finite 
size of the energy source) , within the processing 
performed by the processor 2, and within the display 3; 
limitations arise inherently from, for example, the data 
sampling frequency and display density used, and also 
from the filter-function of the software involved. More 
particularly, and referring to Figure 2, an infinitely 
fine delta function 20 in the object domain is not 
faithfully reproduced in the image domain. Instead, the 
transfer as represented by the arrow 21 results in a 
point-spread function (PSF) or smeared-out image 22 in 
the image domain. Similarly, an infinitely sharp edge- 
response function (ERF) or step 2 3 in the object domain 
becomes through the transfer represented by arrow 24, a 
smeared-out transition 25 of spread represented by 
dimension arrows 26, in the image domain. When two image 
ERFs are close to one another, the smeared-out effects 
run into each other. The consequent deterioration of the 
spatial resolution is often monitored .by the percentage 
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modulation transfer which is given by the ratio, 
expressed as a percentage, of the amplitude of the 
modulation in the image domain to that in the object 
domain. 

The smearing effect becomes more intense as adjacent ERFs 
of discontinuities or contrast profiles get closer to 
each other (or as the spatial frequency of the modulation 
becomes higher) ; this also causes loss of profile height. 
The inherent loss of the spatial resolution (that is, the 
part that is indicated by the smear ed-out effect on the 
corner edge of the ERF) cannot be restored or partially 
restored even by re-scanning the image with an ultra high 
resolution digital scanner system. 

The further processing 5 is operative in accordance with 
the invention to provide accurate image-edge definition 
and location, and to improve spatial resolution in the 
imaging. More especially, in the context of Figure 2, 
the edge position corresponding to the discontinuity or 
step 23 of the object ERF is pin-pointed in the image 
domain from the mid-point of the full-width half-maximum 
(FWHM) of the image PSF; the pin-pointing is to sub-pixel 
accuracy for the image ERF. Low-contrast and A area' 
filtering are used to remove * spurious 7 edges, and sub- 
pixel sampling to detect detail to the resolution of the 
single-pixel modulation. The discontinuity or step 23 of 
the ERF is then restored within the image domain by 
removing the sub-pixel values from outside the optimum 
edge position to compensate for those within. It is to 
be noted that the sub-pixels then become pixels in 
display, and that the enhancement is equivalent to the 
performance of an extra high resolution image transfer 
system . 

As represented in Figure 2 by the arrows 27 and 28, the 
image ERF 2 5 of an infinitely-sharp step 2 3 can be 
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produced by convolution of the image PSF 22 with the 
object ERF 23. In accordance with the present invention, 
de-convolution of the image ERF 25 using sub-pixel 
sampling, represented by the arrow 29, reproduces the 
image PSF 22 in a de-convolution space as image PSF 22 
The image ERF 25 is superimposed on the PSF 22' within 
this space as image ERF 25', and the optimum edge- 
position 3 0 is derived from the mid-point of the FWHM of 
the image PSF 22', and is pin-pointed to sub-pixel 
accuracy . 

For one-dimensional cases, the operation in accordance 
with the invention is relatively simple, as only either 
the x- or the y-profile, that is to say a line spread 
function LSF is involved. But for two-dimensional 
operations, both the x- and y-profiles, and if necessary, 
the xy-diagonal profiles to eliminate any possible 
streakings in the image, may be used; in this case, a 
proper weighting scheme will be required to re-construct 
the image . 

Once the original sharp-edge feature represented by the 
object ERF 2 3 is pin-pointed at the position 3 0 within 
further processing 5, that feature may be restored by 
additional re-processing 6 (Figure 1) , In re-processing 
6, the sub-pixel values occurring x outside' the optimum 
edge-position 3 0 are transferred to compensate those 
v within'. This is illustrated in Figure 2 by arrow 31 
transferring sub-pixel blocks 3 2 from after point 3 0 in 
the image ERF 25', to before it. The re-construction of 
image ERF 25' into image ERF 33 conforming closely in 
configuration to object ERF 2 3 is represented by arrow 
34* Image ERF 33 is displayed in enlarged form in 
display 7 (Figure 1) . 

These techniques enable substantial recovery of the loss 
of spatial resolution in the imaging , without the 
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trade-off loss of other properties such as image noise. 
Furthermore, the enhancement of spatial resolution in 
display 7 reproduces the region of interest selected from 
display 3, without blurring- {or step) effects at the 
profile edge. 

Figure 3 is illustrative of some of the low-contrast 
results provided in practice from MR scanning of a pig's 
brain in fluid. Curve A is the image ERF produced, 
whereas curve B is the line spread function (LSF) 
resulting from de-convolution of curve A carried out in 
processing 5. The optimum edge-position is established 
from the mid-point C of the FWHM of curve B, and the 
additional re-processing 6 is operative by means of sub- 
pixel transfer, to re-construct curve A to conform 
substantially to the edge-feature from which it 
originated in the display 7. 

It is to be noted that whereas curve A is stepped , curve 
B is nonetheless smooth and that mid-point C is located 
to sub-pixel accuracy. Furthermore curve B indicates a 
sensitivity of more than 8:1 between the profile-height 
and background noise. 

Referring again to Figure 2, the complete profile 35 of 
an image within the selected area of interest of display 
3, is built up as indicated by arrow 36 f from the edge- 
position data derived within processing 5. This data 
identifies the location of the point 30, together with 
the locations of all corresponding points derived from 
sampling the multiple x- or y-profiles of the selected 
area of interest. The build up and display of these 
points from the data takes place in display 8 so that a 
substantially true contour 37 for the profile 35 is 
defined. The sharpness of the true contour 37 is in 
contrast to the smeared contour that without de- 
convolution would have been obtained by virtue of the 
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spread (represented by the arrows 2 6) of the relevant 
image ERFs 25. 

A small portion of the contour 37 is shown enlarged in 
Figure 4 and is defined as a best-fit line between 
optimum edge-positions derived respectively from x- and 
y-profiles; the x-profile positions are indicated by 
black dots and the y-profile positions by white dots. 
The closeness to one another of corresponding positions 
in the x- and y-profiles is indicative of the accuracy to 
sub-pixel level achieved. The smear-out that would have 
been manifest in the image-profile contour if the de- 
convolution technique were not used, would have extended 
throughout the space bounded by the dashed lines 3 8 and 
39; these boundaries are indicated in dashed line in the 
representation of profile 35 in Figure 2. 

The accurate definition of the image contour 37 derived 
in the display 8 allows correspondingly accurate 
determination in calculation 9 of the area within that 
contour; the volume involved can also be derived from 
successive slices scanned. The determination of area and 
volume is especially useful for diagnostic and accurate 
assessment 10 of the size of a tumour or lesion before 
and after treatment. It is similarly useful for 
assessment of arterial dimensions in angiography. 

Moreover, the accurate definition of the image contour 37 
derived in the display 8, is particularly useful for 
segmenting anatomical structures for diagnostic and 
treatment planning 11. Furthermore, the ratio of 
intensities of two scans are derived by processing 12 to 
derive values of relaxation times T x and T 2/ and of proton 
density. The values are then represented in display 13 
within the boundary of the image contour, utilising 
standardisation data derived from a couch-top 14 used 
within the scanner 1. The couch-top J.4, which also 
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provides landmarks for determining position coordinates, 
has the form shown in Figures 5 and 6, and will now be 
described. 

5 Referring to Figures 5 and 6, the couch-top 40, which is 
of polystyrene bubble foam, has the form of a flat slab, 
and is supported on a curved foam base 41 that fits to 
the scanner table of the MR installation. Two zig-zag 
runs of tubing 42 are embedded within the top 40 to 
10 extend along the length of the couch either side of 
centrally-embedded straight tubing 43. 

The tubing 42 of each zig-zag run is of double-bore 
rectangular cross-section, whereas the tubing 4 3 is of 

15 single-bore cross-section. The five bores defined by the 
array of tubing 42 and 43 may be filled respectively with 
the five MR solutions S 0 to S 4 of Table I, for 
standardisation and calibration purposes. The four 
MnCl 2 .4H 2 0 solutions, S } to S 4/ cover the full range of 

20 values of T x and T 2 for anatomical tissues, and the fifth 
solution, S 0f of CuS0 4 .5H 2 O, is nominally equivalent to 
"loosely bound water". 



Table I 



Solution 




T, 

at 0,5T 


T 2 

at 0.5T 


So 


1.25 g/1 CuS0 4 .5H 2 O 


200 ms 


200 ms 


s, 


3.41X10 16 Mn +2 ions/ml 


840 ms 


300 ms 


s 2 


1.15xl0 17 Mn +2 ions/ml 


440 ms 


120 ms 


s 3 


2.30X10 17 Mn +2 ions/ml 


2 50 ms 


60 ms 


s 4 


4-37X10 17 Mn +2 ions/ml 


150 ms 


3 0 ms 
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The tissue types revealed by the T,, T 2 and proton density 
values in display 13 are determined by processing 15 from 
look-up tables, and tissue densities are assigned within 
the image-contour boundaries in display 16. The images 
of display 16 are furthermore corrected empirically for 
geometry distortion in accordance with data supplied from 
memory 17. The data stored in memory 17 is derived using 
a set of drum phantoms of the spider-web form, and 
correction for geometry distortion is realistically 
effective principally because of the accuracy and spatial 
resolution with which image-contours are defined. 

The tissue types assigned to the corrected images are 
utilised in display 16 through the Bulk Heterogeneity 
Correction method described by Richard A, Geise et al , 
Radiology, 124:133-141, July, 1977, to establish for each 
image a normalised tissue density value; the up-datable 
look-up table for this is stored in memory 18. 
Accordingly, the display 16 when used at step 19 in 
conjunction with the positional datums derived from the 
couch 14, has all the tissue contours accurately mapped 
out with their respective tissue densities and locations. 
This establishes an accurate and readily-usable, stand- 
alone basis for diagnosis and treatment planning, and 
enables a true three-dimensional assessment and plan to 
be made when both orthogonal and oblique MR images are 
involved. 

Although the method and system of the invention have been 
described above in the medical context they are 
applicable more widely than this, for example, in 
engineering, in physical science and in the field of 
instrumentation generally. Moreover, the method and 
system is not limited to MR imaging, but may be utilised 
where other forms of imaging are involved. The steps and 
structure represented in Figure 1 by x boxes' 1 to 11 are 
just as applicable to computer assisted tomography (CT) , 
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as they are to MR imaging. Other forms of imaging to 
which the invention is applicable include X-ray 
radiography, film- or print-image transformation to 
digital form, digital X-ray f luorography , ultra-sound 
5 imaging, nuclear medicine, positron emission tomography 
(PET) and other camera or imaging. The technique is 
particularly suitable for use in X-ray digital 
f luorography , in which small structures under study are 
highlighted by injection of contrast liguids; the small 
10 structures may also be isolated from surrounding 
interfering effects by using an image subtraction 
technique. 

The inherent resolutions of X-ray radiography, 
15 ultra-sound imaging, nuclear medicine, and PET scanning 
are relatively low, and some are used for real-time 
study. Only the individual still frame or hard-copy 
images may be re-processed. 

20 In the context of engineering, physical science and the 
field of instrumentation, the invention is applicable to 
one-dimensional imaging as used, for example, in regard 
to bar-code patterns, the spectrum of DNA analysis, iris 
patterns of eyes (for example, for identification 

25 purposes in commercial banking) , finger-print 

identification, and emission spectroscopy. The invention 
is also applicable to two-dimensional imaging, for 
example, in relation to images obtained by satellite or 
pattern recognition, or from a surveillance camera or 

3 0 during laboratory experimentation. As a general matter, 
the invention is applicable where there needs to be 
accurate determination of the edge position in an image 
versus the true object-edge position, for the purpose, 
for example, of measurement of the positional 

35 displacement between object and image, distortion 
correction and manufacturing control. 
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As a further example of application of the present 
invention, a method and system that uses CT and MR 
imaging in conjunction with one another, will now be 
described. 

The major contribution to the magnetic-resonance (MR) 
signal comes from the abundant protons content of water 
molecules and protein. It is a quantum process at the 
Larmor frequency according to the magnetic field in use. 
The ^-weighted' and ^T 2 -weighted' MR signals from 
protons provide contrast numbers that are relative in 
scale, whereas in CT, the X-ray absorption is a 
polychromatic attenuation process affected by the 
electron densities of all the atoms presented within the 
X-ray beam. There is no equation to correlate the CT 
number (or the linear attenuation coefficient, electron 
density, or tissue density) with the MR-contrast numbers; 
no direct calibration between the two types of signal is 
possible. This lack of correlation is confirmed by 
consideration of bone and air which are at opposite ends 
of the CT contrast (absolute) scale using water as the 
base-line reference, but which are at the same end of the 
MR-image contrast (relative) scale owing to their common 
low proton-population . 

The lack of correlation between the CT and MR signals 
acts against their use in combination for imaging 
purposes, but the present invention provides a method and 
system by which the advantages of each may be utilised to 
improve image resolution and contrast information. 

In the latter regard, CT provides a high spatial 
resolution but only in regard to view normal to the 
transverse slice. Resolution for all re-constructed 
non-transverse planes is poor owing to the need to use 
elongate voxels to improve signal-to-noise ratio. Also, 
the partial-volume effect of using elongate voxels may 
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give rise to detection errors at the thin edge of a 
contrast profile of a lesion. MR, on the other hand, can 
give the same high degree of spatial resolution viewed in 
the normal direction to any image-slice plane, and can 
5 also provide isotropic resolution with cubic-voxel volume 
imaging. 

To this end, multiple-slice transverse CT scans are 
collected across a section of the volume of interest in a 

10 patient or other subject. Corresponding multiple-slice 
transverse MR scans of the same volume are also 
collected. The slice thickness of the latter scans may 
be one half, or smaller, of the thickness of the CT 
slices, and may be collected two~dimensionally or three- 

15 dimensionally . The patient or other object scanned is 

constrained throughout on a couch that provides landmarks 
with respect to a coordinate reference arrangement on the 
couch-top. This is to ensure the reproducibility of, for 
example, anatomical positions and features to the first 

20 order accuracy for the corresponding CT and MR scans, and 
possibly for radiation treatment to be made. 

The respective transverse planes of the CT and MR images 
are processed individually in the method and are matched 

25 with, one another in a de-convoluted space for the CT and 
MR images. The two sets of de-convoluted maps are then 
merged together to a second order of accuracy in order 
that the CT numbers may be transferred over to replace 
the corresponding MR contrast numbers. Once this has 

3 0 been achieved, non-transverse (or oblique) planes can be 
obtained from the two-dimensional MR images, or from the 
re-arrangement of the corresponding voxels of the three- 
dimensional volume images; where two-dimensional MR is 
used, a further step of contrast transformation may be 

3 5 required. 
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The respective transverse planes of CT and MR images are 
processed individually by using x boundary' or 
x finger-print' matching techniques in a de-convolution 
space for the CT & MR images. In these transverse image 
5 planes, in particular when used medically, the 

skin-contour features along the sides of the patient may 
be best used for second-order alignment and matching 
purposes, as they are less affected by patient-movements* 
Transitional error may be readily corrected with respect 

10 to the coordinate positions of a rectangular tubing 
system embedded in the couch-top (for example that 
described above with reference to Figures 5 and 6) . The 
processed data may then be used for a % diagnostic and 
statistics software package' of CT image versus MR image 

15 for their exactly corresponding transverse slice(s) / and 
an associated x statistical package 7 for accurate 
computation of the x true' area, and then the 'true' 
volume, of a lesion or tissue profile or contour. 

20 The two sets of de-convoluted maps may also be merged 

together to a second order of accuracy in order that the 
CT numbers may be transferred over to replace the 
corresponding MR contrast numbers. Once this has been 
done, non-transverse (or oblique) planes are obtained 

25 from the two-dimensional image or from the re-arrangement 
of the corresponding voxels of the three-dimensional 
volume; in the two-dimensional MR case, a further step of 
contrast transformation is required. The transferred 
contrast data may then be used in a three-dimensional 

30 radiotherapy treatment planning software package for an 

in-plane, oblique-image pseudo-three-dimensional approach 
using MR images. 

Software required for de-convolution processing of image 
35 data according to the invention may be implemented in 

conjunction with a least-squares curve-fitting method. A 
mathematical model of the method, frop which the required 
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software can be readily developed, will now be given in 
relation to Figure 7 which indicates the convolution of 
an LSF L(x) with a step function SF having a step-down 
edge SDE, to produce an image ERF E(x) . Considering the 
values of E(x) resulting from the convolution at points 
x=l to x=n of SF: 



at x=l; 



10 at x=2; 



£l,<x) Ax = E(x) 1 



J^L(x)Ax - Y^L(x)Ax - E(x) 



at x=3 



J^L(x) Ax - Y^L(x) Ax * E(x) 



a-2 



until at x=n: 



£ L (x) Ax - £ L (x) Ax = E(x) n 



For de-convolution, that is to say, for the ERF E(x) to 
manifest the LSF L(x) : 

a 

E{X) x - E{X) 2 = L Ak = L{x) (a-l)-a Ax 

a-I 

{E(x) x - E(x) 2 )/Ax = L(x) (a . 1)-a 

thus: 

iE(x) 2 - E(x)J/Ax = Z,U> (a - 2) - (a . 1} 

leading to 



20 
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Thus, the same shape of LSF is recovered by the 
de-convolution process independently of the sense of the 
ERF; it does not matter whether a roll-up or a roll-down 
ERF of a positive or a negative contrast profile, is 
5 involved. This is an important property as a positive 
contrast contour will have roll-up (from low to high) 
ERFs at both ends whereas a negative contrast contour 
will have roll-down (from high to low) ERFs at both ends. 

10 In the practical implementation, the LSF is derived from 
the ERF by de-convolution using a running filter. This 
enhances the accuracy of the method in overcoming the 
problem of noise that affects the digitised pixel values 
of the image. Use is made of a least-squares fitting of 

15 a set of data points that advances along the whole length 
of the function from one point to another. 



20 



Assuming that y = a 0 + a,x + a 2 x 2 represents the ERF curve 
a five-point or seven-point fit is used, and the normal 
equation becomes : 



Ei Hx i Ex| 
Hx i Ex? Ex. 3 
Ex] Ex, 3 Ex! 



aO 
al 
a2\ 



Ex i y i 
Ex? yi 



25 



where: i is 1+n, ... 5+n; ... until m-5, ... m 
or 

i is l+n, ... 7+n; ... until m-7 , ... m 
n is 0 , 1 , 2 , ... ; and 

1, ... m is the span of the ERF profile. 



The solution may be derived from either: 



aO 


El 


Ex 


Ex 2 


-1 


Ey- 


al -- 


= Ex 


Ex 2 


Ex 3 




Exy 


a2 


Ex 2 


Ex 3 


Ex 4 




Ex 2 y 
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aO 

Ey Ex Ex 2 
Exy Ex 2 Ex 3 

Ex 2 y Ex 3 Ex 4 



al 

El Ey Ex 2 
Ex Exy Ex 3 

Ex 2 Ex 2 y Ex 4 
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a2 

" El Ex Ey 
Ex Ex 2 Exy 
Ex 2 Ex 3 Ex 2 y 



1 

El Ex Ex 2 
Ex Ex 2 Ex 3 

Ex 2 Ex 3 Ex 4 



For both of these equations to be valid 



Ex 



Ei E 
Ex Ex 2 Ex 3 
|Ex 2 Ex 3 Ex 4 



* o 



5 The gradient at dy/dx (3+n) or dy/dx (4+n) can then be derived 
and plotted against x for the LSF profile. 



The graph of dy/dx against x gives the LSF profile. The 
peak of this profile is located centrally of the mid- 
10 points of the ascending and descending limbs of the 

graph. These points define the extremes of the full- 
width-half -maximum (FWHM) of the profile and the mid- 
point of this is determined with an accuracy of sub-pixel 
level owing to the ^average' effect. 

15 

The point spread function (PSF) is the two dimensional 
profile which may be derived, in practice, from the two 
corresponding LSFs orthogonal to one another within an 
image plane. The peak position of the PSF profile is, 
20 therefore, from the x mean' or x cross-over' of the two 
peaks or the two LSF profiles. The PSF is obtained in 
practice from two orthogonal axes in a two-dimensional 
plane* 



25 The generation of the LSF (or PSF) is, after phase- 
reversal correction, independent of the roll-up or 
roll-down nature of ERFs at the edges of the contrast 
contour. In other words, it is independent of the sense 
and the absolute value of the contrast' numbers within the 
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contour. The peak position of the LSF (or PSF) is the 
central half-way (50%) point of the roll-up or roll-down 
ERF which is the optimum position for true-edge 
definition. 



in 
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Claims: 

1. A method of imaging in which a de-convolution 
process is applied to the image-domain results of an 
object-scan to derive therefrom the respective point- or 
line-spread function of one or more object- 
discontinuities, and to derive from said function the 
location in the image domain of the respective 
discontinuity, wherein the de-convolution process is 
carried out using sub-pixel sampling* 

2. A method according to Claim 1 wherein the location 
of the respective discontinuity is derived from the mid- 
point of the full-width half -maximum of said function. 

3. A method according to Claim 1 or Claim 2 wherein 
said function is correlated with the image-domain results 
of said transfer for enhancement of spatial resolution of 
the imaging of the one or more discontinuities. 

4. A method according to Claim 3 wherein the 
enhancement of spatial resolution of the imaging of the 
one or more discontinuities involves transfer of sub- 
pixels within the image-domain results of the respective 
one or more discontinuities, the sub-pixels being 
transferred within their respective image-domain results 
from one side to the other of said location for edge- 
image definition. 

5. A method according to any one of Claims 1 to 4 
wherein the de-convolution process is carried out using 
least-squares running filtering. 

6. A method according to any one of Claims 1 to 5 
wherein an edge-contour of the object is defined in the 
image domain using de-convolution processing as 
aforesaid. 
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7. A method according to Claim 6 wherein the area 
and/or volume of the object-image within the edge-contour 
is determined. 

8. A method according to Claim 6 or Claim 7 wherein the 
object-scan is a magnetic resonance (MR) scan, values of 
relaxation times T x and T 2 are derived for the object- 
image within said contour , and these values are used to 
identify from stored data, types of tissue or other 
material involved in the scanned object, 

9. A method according to Claim 8 wherein density values 
for the identified tissue or other material types are 
derived from further stored data. 

10. A method of imaging according to any one of Claims 1 
to 9 wherein corresponding computed tomography (CT) and 
magnetic resonance (MR) scans of the same part of an 
object are derived, the scans are related to one another 
for correlation of one to the other positionally with 
respect to said part using the de-convolution process, 
and imaging of said part of the object is provided in 
accordance with the MR scan as modified spatially in 
dependence upon the CT contrast numbers applicable to the 
corresponding, correlated positions of the CT scan. 

11. A method according to any one of Claims 8 to 10 
wherein geometric correction is applied to the imaging 
derived from the MR scan, in accordance with stored data. 

12. An imaging system comprising means for performing a 
de-convolution process on the image-domain results of an 
object-scan to derive therefrom the respective point- or 
line-spread function of one or more object- 
discontinuities , and means to derive from said function 
the location in the image domain of the respective 
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discontinuity, wherein the de-convolution process is 
carried out using sub-pixel sampling. 

13 . A system according to Claim 12 wherein the location 
of the respective discontinuity is derived from the mid- 
point of the full-width half -maximum of said function. 

14. A system according to Claim 12 or Claim 13 wherein 
said function is correlated with the image-domain results 
of said transfer for enhancement of spatial resolution of 
the imaging of the one or more discontinuities. 

15. A system according to Claim 14 wherein the 
enhancement of spatial resolution of the imaging of the 
one or more discontinuities involves transfer of sub- 
pixels within the image-domain results of the respective 
one or more discontinuities, the sub-pixels being 
transferred within their respective image-domain results 
from one side to the other of said location for edge- 
image definition. 

16. A system according to any one of Claims 12 to 15 
wherein the de-convolution process is carried out using 
least-squares running filtering. 
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IMAGING 
ABSTRACT OF THE DISCLOSURE 



A de-convolution process is applied to an MR, CT or other image (25) of a 
scanned-object (23) to derive the point-spread function (22') at an object-edge and 
to pin-point from the mid-point of its full-width-half-maximum FWHM, the location (30) 
of the true image-edge. With the object-image (25') overlying the PSF function (22') 
in the de-convolution space, sub-pixels which follow location (30) are transferred to 
before it to re-construct the image-edge (25') for sharper conformity to the object- 
edge (23). Sharp definition of image-contour (37) facilitates accurate determination 
of area and volume of image profiles (35) and their segmentation. The accurate 
image-edge definition enables viable correction of geometrical distortion in stand- 
alone MR diagnosis and treatment planning. 
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application and transact all business in the Patent and Trademark Office connected therewith, (list 
name(s) and registration number(s)) 

Anthony G. M. Davis Registration No. 27,868 

Michael J. Bujold Registration Nd~323^ 

Scott A. Daniels Registration No?42^462 

□ Attached as part of this Declaration and Power of Attorney is the authorization of the above- 
named attorney(s) to accept and follow instructions from my representative(s). 




Send Correspondence to: 

^ayi^^ujoldJ , . L. L. C. 
Fo urth Floor 

500 N. Commercial Street 
Manchester, NH 03101-1151" 



Direct Telephone Calls to: 
(603) 624-9220 



Direct Telefaxes to: 
(603) 624-9229 



ACKNOWLEDGMENT OF REVIEW OF PAPERS AND DUTY OF CANDOR 



I/We hereby state that l/we have reviewed and understand the contents of the above identified 
specification, including the claims, as amended by any amendment referred to above. 

I/We acknowledge the duty to disclose to the United States Patent Office all information which 
is known to be material to patentability of this application as defined in § 1 .56 of Title 37 of the Code of 
Federal Regulations. 

PRIORITY CLAIM 

[/We hereby claim foreign priority benefits under Title 35, United States Code, § 1 1 9 of any foreign 
application(s) for patent or inventor's certificate or of any PCT international application(s) designating at 
least one country other than the United States of America listed below and have also identified below any 
foreign application (s) for patent or inventor's certificate or any PCT international application(s) designating 
at least one country other than the United States of America filed by me/us on the same subject matter 
having a filing date before that of the application(s) of which priority is claimed. 



EARLIEST FOREIGN APPLICATION(S), IF ANY FILED WITHIN 12 MONTHS 
(6 MONTHS FOR DESIGN) PRIOR TO THIS U.S. APPLICATION 



COUNTRY 


APPLICATION NO. 


DATE OF FILING 
(day, month, year) 


PRIORITY CLAIMED 
UNDER 37 USC 119 


Great Britain 


9822397.7 


October 15, 1998 


■YES DNO 


Great Britain 


9825165.5 


November 18, 1998 


■YES UNO 


Great Britain 


9902332.7 


February 2, 1999 


■YES DNO 








□YES DNO 








□YES DNO 



ALL FOREIGN APPLICATION(S), IF ANY FILED MORE THAN 12 MONTHS 
(6 MONTHS FOR DESIGN) PRIOR TO THIS U.S. APPLICATION 



q I/We hereby claim the benefit, under 35 U.S.C. 119(e), of any United States provisional 
application(s) listed below. 



Application Number(s) 


Filinq Date fMM/DD/YY) 


□ Additional provisional 
application numbers are listed 
on a supplemental priority data 
sheet PTO/SB/02B attached 
hereto. 







DECLARATION 

I/We hereby declare that all statements made herein of my/our own knowledge are true and that 
all statements made on information and belief are believed to be true; and further that these statements 
were made with the knowledge that willful false statements and the like so made are punishable by fine 
or imprisonment, or both, under Section 1001 of Title 18 of the United States Code, and that such willful 
false statements may jeopardiz^he validity of the application or any patent issued thereon. 

Full name of sole inventor ! Kui Ming CHUI _ >a a 

k. ~. ~ An*-, 2* ei 



Inventor's signature:. 



Residence: 8 Gilbev Close. Ickenham. Uxbridge, Middlexes UB10 8TD. Great Britain 

Post Office Address: Same as above Country of Citizenship: Great Britain 



